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ABSTRACT 

■ The early evolution of the magnetic field and angular momentum of newly 

^ ■ formed protostars are studied using three-dimensional resistive MHD nested grid 

CN ■ simulations. Starting with a Bonnor-Ebert isothermal cloud rotating in a uniform 

magnetic field, we calculate the cloud evolution from the molecular cloud core 
[ric — lO^cm"^, r = 4.6 x 10^ AU, where and r are central density and radius, 
respectively) to the stellar core {ric — lO^^cm"^, r ~ l-R©). The magnetic 
field strengths at the center of clouds with the same initial angular momentum 
^ . but different magnetic field strengths converge to a certain value as the clouds 

^ . collapse for ric ^ lO^^cm"^. For 10^^ ^ '^c ^ lO^^cm"^, Ohmic dissipation 

largely removes the magnetic field from a collapsing cloud core, and the magnetic 
^ . field lines, which are strongly twisted for ric ^ lO^^cm"^, are de-collimated. 

. The magnetic field lines are twisted and amplified again for ric > lO^^cm"^, 

because the magnetic field is recoupled with warm gas. Finally, protostars at 
'"q_|. their formation epoch {ric — 10^^ cm'^) have magnetic fields of ~0.1-1 kG, which 

is comparable to observations. The magnetic field strength of a protostar depends 
slightly on the angular momentum of the host cloud. A protostar formed from 
d '. a slowly rotating cloud core has a stronger magnetic field. The evolution of the 

angular momentum is closely related to the evolution of the magnetic field. The 
^ ! angular momentum in a collapsing cloud is removed by magnetic effects such as 

magnetic braking, outflow and jets. The formed protostars have rotation periods 
of 0.1-2 days at their formation epoch, which is slightly shorter than observations. 
This indicates that a further removal mechanism for the angular momentum, such 
as interactions between the protostar and the disk, wind, or jets, is important in 
the further evolution of protostars. 

Subject headings: ISM: clouds: ISM: magnetic flelds — MHD — stars: formation — 
stars: rotation 
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Introduction 



The Lorentz and centrifugal forces play important roles in the star formation process. 
While gravity and thermal pressure are isotropic forces, the Lorentz and centrifugal forces 
are anisotropic forces and are closely related to disk formation, outflow, and jets in col- 
l apsing clouds. Molecu l ar clouds hav e ~ 2% rotational energy against gravitational energy 



( iGoodman et al. Il993l: ICaseUi 



20021 ). while the magnetic energy is comparable to the grav- 



itational energy (jCrutcherl Il999l ). By considering the conservation of magnetic flux and 
angular velocity, it is seen that the rotation and magnetic field in a cloud gradually increase 
as the cloud collapses. However, the magnetic field strength and angular velocity of observed 
protostars indicate that neither the magnetic flux nor angular momentum is conserved in 
collapsing clouds. In general, these anomalies are called the 'magnetic flux problem' and 
'angular momentum problem.' The former problem refers to the fact that the magnetic flux 
of a molecular cloud is much larger than that of a protostar with equivalent mass. The latter 
problem is that the specific angular momentum of a molecular cloud is much larger than that 
of a protostar. These problems imply that there must be mechanisms removing magnetic 
flux and angular momentum from a cloud core. In a collapsing cloud, these two problems 
are mutually related. Namely, the angular momentum is removed by magnetic effects (i.e., 
magnetic braking, outflow, and jets), while the magnetic field is amplified by the shearing 
motion caused by cloud rotation. Hence, the magnetic field and rotation cannot be treated 
independently in considering the magnetic flux and angular momentum problems. In addi- 
tion, it is difficult to treat the evolution of the magnetic field and rotation analytically since 
the density {n ~ lO^cm"^) and scale {R ~ 10^ AU) of molecular clouds are very different 
from those of protostars (n ~ lO^'^cm"^, R ~ l-R©)- Therefore, numerical simulation is 
needed to study the magnetic field and rotation of protostars formed from molecular clouds. 

Angular momentum is removed from a collapsing c l oud b y mag netic braking and outflow, 
as shown bv iBasu fc Mouschoviad (119941 ). iTomisakal (120021 ) . and iMachida et al. I (l2005al ). 
Machida et al. I (j2005al ) found that about 70% of the total angular mome ntum is removed 
from the collapsing cloud core in the isothermal phase {n < lO^^cm"^). ITomisakal (120021 ) 
shows that ~99% of the total angular momentum is transferred by outflow in the adiabatic 
phase (lO^^cm^^ ^ ""^ ^ lO^^cm"^). Thus, previous studies show that angular momentum 
is effectively transferred by magnetic effects. In contrast, there are relatively few studies of 
the evolution and removal pr ocess of magnet i c flux in collapsing cloud cores, especially for 
high density [n > 10^^ cm~^). iNakano et al. I (|2002| ) analytically investigated the dissipation 
process of the magnetic field in a collapsing cloud core and showed that it is expected that 
the magnetic fiux is removed from the cloud core largely for 10^^ cm~^ ^ ^ 10^^ cm~^. 



To investigate the evolution of the magnetic field and rotation (or the magnetic fiux 
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and angular momentum problems), it is necessary to calculate the cloud evolution of a 
cloud from the molecular cloud core {n ~ lO^cm"^) to protostar formation [n ~ 10^'* cm^^) 
while accounting for the magnetic field and rotation. The evolution of a magnetized cloud 
up to a centra l dens it y of n 10^^ cm 



has been i nvest igated by Hosking &: WhitworthI 

fl2004[).lzieglCTl fl2005yMatsumoto fc Tomisaka booj ). and lMachida et al. I J2004j . boOSal lbh. 
Tomisakal (|2002( ) and lBanerjee fc Pudrita (120061 ) calculated cloud evolution up till the forma- 
tion of a protostar {n ~ 10^^ cm~^), and showed that the magnetic field plays important roles 
in the star- formation process. However, they adopted an ideal MHD approximation, which is 
valid in low-density gas regions [n < 10^^ cm^^), but is not valid in high-density gas regions 
{n > lO^^cm"'^). A significant magneti c flux loss occurs for lO^^cm"^ ^ ""^ ^ lO^^cm^^ by 
Ohmic dissipation (INakano et al. II2002I ) and hence these studies overestimate the magnetic 
flux in a collapsing cloud, especially in high-density gas regions. 



Following iMachida et al. I (j2006bl ). we report detailed results of three-dimensional non- 
ideal MHD nested grid simulations. In this paper, we calculate the cloud evolution from 



the cloud core 



.rir. 



10^ cm 



4.6 X 10^ AU) to protostar formation {ric 



10^2 cm-\ 



rc — 1 Rq) and discuss the magnetic field and angular momentum of the formed protostar. 
The structure of the paper is as follows. The frameworks of our models are given in §2 
and the numerical method of our computations is shown in §3. The numerical results are 
presented in §4. We discuss the magnetic field and angular momentum of a protostar and 
compare our results with observations in §5. 



2. Model 

2.1. Basic Equations 

To study the cloud evolution, we use the three-dimensional resistive MHD nested grid 
code. We solve the resistive MHD equations, including self-gravity: 

f^+V-{pv) = 0, (1) 

dv 1 
p— + p{v ■V)v = -VP - —B X {V X B) - pV(j), (2) 
ot Air 

f) Ft 

— = S/x{vxB)+ 7]V'^B, (3) 
VV = AttGp, (4) 



where p, v, P, B, rj, and denote the density, velocity, pressure, magnetic flux density, 
resistivity, and gravitational potential, respectively. The last term in equation (3) denotes the 
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Ohmic dissipation. Although the dissipation term is expressed exactly by —V x (r^V x B), we 
simplify it as shown in the equation (3) . This simplification allows us to reduce computational 
costs considerably, and we can follow the collapse of the cloud up to the stages of the stellar 
core formation. The deviation due to this simplification is examined by test calculations of 
typical models; this simplification does not affect evolution of the cloud considerably, and 
the magnetic field strength deviates by a factor of two at most at the central region. 



To mimic the temperature evolution calculated by iMasunaga fc Inutsuka I (120001 ) , we 
adopt the piece-wise polytropic equation of state: 



P 



7/5 

7/5 



Pc 

Pd\' I JL 
Pc J \Pd 

7/5 / \ 1.1 



P < pc, 
Pc< P < Pd, 

Pd < P < pe, 



(5) 



5/3 



P > Pe, 



190 ms \ p, 



3.84 X lO^^^g cm-3 



10" cm~^) 



pd 



3.84 X 10~*^g cm-3 



where Csfi = 

{ud = lO^^cm"^), and pe = 3.84 x lO^'^g cm~^ (rig = lO^^cm"'^). For convenience, we 
define 'the protostar formation epoch' as that at which the central density reaches ric = 
We also call the period for which ric < 10^^ cm~^ 'the isothermal phase', the 



10^1 cm-3 



< ric < 10^^ cm 'the adiabatic phase', the period for which 



period for which 10^^ cm" 
10^^ cm^'^ < ric < 10^^ cm~'^ 'the second collapse phase', and the period for which ric > 
10^^ cm~^ 'the protostellar phase.' We, therefore, ignored the effect from heating by Ohmic 
resistivity, because we adopted a simple equation of state. Plasma beta in the collapsing 
clouds maintains (3p > 1 (see §4). Thus, even if all the magnetic field energy is converted 
into the thermal energy, the gas temperature increases by a factor of two. However, since 
plasma beta is /3p > 10 in any model when Ohmic dissipation becomes effective (see §4.1), 
we can safely neglect the heating effect by Ohmic resistivity. 



In this p aper, we adopt a spherical cloud with a critical Bonnor-Ebert (lEbertl Il955 



Bonnorlll956l ) density profile pbe as the initial condition. Initially, the cloud rotates rigidly 
with angular velocity Qq around the 2;-axis and has a uniform magnetic field Bq parallel to 
the z-axis (or rotation axis). To promote contraction, we increase the density by a factor / 
(density enhancement factor), as 



p(r) 



pBEir) f for r < 
Pbe(^c) / for r > 



(6) 



where r and Rc denote the radius and the critical radius for a Bonnor-Ebert sphere, re- 
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spectively. We adopt density enhancement factors of / = 1.2 and 1.4 (see Table [T]) B We 
assume Pbe(O) = / x 3.841 x 10~^°g cm~^, which corresponds to a central number den- 
sity of ricfi = / X lO^cm"^. Thus, the critical radius of a Bonnor-Ebert sphere Rc = 
6.45c,,o[47rGpB£(0)]-^/2 corresponds to = 4.58 x lO^AU for our settings. The initial 
model is characterized by three nondimensional parameters: a, u, and c^. The magnetic 
field strength and rotation rate are scaled using a central density po = Pbe(O)/ as 

« = 5oV(4vrpo<o), (7) 

u = no/{^7rGpoy/'. (8) 

The parameter represents the degree of resistivity (for details, see the next section). We 
calculated 33 models by combining these parameters, which are listed in Table [H The model 
parameters a, uj, and Cr,; density enhancement factor /; magnetic field Bq; angular velocity 
Oq; total mass M inside the critical radius (r < Rc); and the ratio of the thermal, rotational, 
and magnetic energies to the gravitational energy ao, Po, and 70 are summarized in this table 
I. 



2.2. Resistivity and Magnetic Reynolds Number 



To describe realistic evolution of magnetic field in the protostar formation, we should 
take into account the non-ideal MHD effects of weakly ionized molecular ga s. The detailed 
physi cal processes in the problem at hand are studied by many authors (e.g.. iNakano et al. 
2OO2I . references therein). In general, the ambipolar diffusion is slow but important in the 
low density phase, and Ohmic dissipation dominates in the high-density phase. In the in- 
termediate density phase, t he Hall term effect can be also important depending on the size 
distribution of dust grains (IWardle and Ng 1119991 ). Note, however, that Ohmic dissipation 
is the most efficient mechanism for the dissipation of magnetic field in the magnetically 
supercritical cloud core, and we are mainly interested in this dynamically contracting gas. 
Therefore we model the dissipation of magnetic field only by the effective resis tivity in the in 
duction equation ([3]). We quantitatively estimate the resistivity r} according to lNakano et al. 



^ The density enhancement factor is related to the stabihty of the initial cloud. The cloud is more unstable 
with larger /. Comparativel y stable clouds are con sidered in this paper. However, cloud evolution hardly 
depends on / (for details see lMachida et al. Il2006al ). 

^ Denoting the thermal, rotational, magnetic, and gravitational energies as U, K, M, and W, the relative 
factors against the gravitational energy are defined as ao — U/\W\, f3o — K/\W\, and 70 = Af/|iy|. 
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( 120021 ) and assume that 1] is a. function of density and temperature as 

cm'^s-\ (9) 



_ 740 T 
^ ~ xTV lOK 



1 ' ^ 

tanii 



10^^ cm 



where T and n are the gas temperature and number density, and Xe is the ionization degree 
of the gas which describes as 



Xe = 5.7 xlQ-U^^] . (10) 

We added the second term in the right-hand side of equation to smoothly dechne the 
diffusivity at n ~ 10^^ cm~^, which means Ohmic dissipation becomes ineffective for n > 
10^^ cm~'^. When the gas density reaches n ~ 10^^ cm~^, the temperature reaches T ~ 2000K. 
Thus, the thermal ionization of alkali metals reduces the resistivity and the magnetic field is 
coupled with the gas n > 10^^ cm~^. We are using a barotropic equation of state, and hence, 
the temperature is a function of the density. Therefore, t] can be expressed as a function of 
density, represented by the thick line in Figure [H To take into account the uncertainty of 
the effective resistivity, we parameterize the maximum value of rj, keeping the shape of the 
function, and define c,, as 

V = Cr,Vo{p), (11) 

where ?7o(p) is a function of the central density and corresponds to the thick line in Figure [TJ 
1] corresponds to ?7o(p) when Cr, = 1. In this paper, we adopt = 0, 10"'^, 0.01, 0.1, 1, and 
10, as listed in Tabled! to investigate the dependence on r], and hence, the effect of Ohmic 
dissipation. Since the second term on the right-hand side in Equation Q vanishes in models 
with = 0, these models obey the ideal MHD approximation. In this paper, we call models 
having Crj ^ 'non-ideal MHD models' and models having Crj = 'ideal MHD models.' 

We analytically estimate the magnetic Reynolds number for models with different 
(Rem = VfXjT]'^; thin lines in Fig. [T]), where Vf = [(4/3)7rG'A|pc]^^^ is the free-fall ve- 
locity and Xj = (rrcl/Gpcy^'^ is the Jeans length. Magnetic dissipation is effective for 
2 X 10^^ cm"'^ ^ ""-c ^ 6 X 10^^ cm~^ in mode ls with Crj = 1 (thin solid line), which cor- 
responds to the results of iNakano et al. I (|2002| ). while it is effective for 2 x 10^^ cm~^ ^ ^ ^ 



2 X 10^^ cm~^ in models with = 0.01. In contrast, at density of ~ 10^^ cm~^ for example, 
a mean free path of main particles is about 1cm. Thus, (hydrodynamic) Reynolds number 
is extremely large, and magnetic Prandtl number is correspondingly small. Note that the 
viscous dissipation is important only at the shock front, and we can safely neglect the effect 
of (hydrodynamical) viscosity in the calculation with our shock-capturing numerical scheme. 
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Numerical Method 



We adopt the nested grid method (for details, see iMachida et al. I l2005al l2006al ) to 



obtain high spatial resolution near the center. Each level of a rectangular grid has the same 
number of cells (64 x 64 x 32), although the cell width h{l) depends on the grid level /. The 
cell width is reduced by a factor of 1/2 as the grid level increases by 1 (/ ^ / + We assume 
mirror symmetry with respect to z = 0. The highest level of a grid changes dynamically. We 
begin our calculations with four grid levels (/ = 1, 2, 3, 4). The box size of the initial finest 
grid / = 4 is chosen to be 2Rc, where Rc denotes the radius of the critical Bonnor-Ebert 
sphere. The coarsest grid, / = 1, therefore, has a box size of 2^i?c- A boundary condition 
is imposed at r = 2^ R^, such that the magnetic field and am bient gas rotate at an angular 
velocity of Qq (for details, see iMatsumoto fc Tomisakall2004l ). A new finer grid is generated 
whenever the minimum local Aj becomes smaller than 8 /i(/max)- The maximum level of grids 
is restricted to /max = 30. Since the de nsity is highest in the finest grid, the generation of a 
new grid ensures the Jeans condition of iTruelove et al. I (119971 ) with a margin of sa f ety fa ctor 



of 2. We adopted the hyperbolic divergence B cleaning method of iDedner et al. I (120021 ) 



Results 



Starting from the number density nc = / x 10^ cm ^ (/ = 1.2 or 1.4), we calculate the 

lO^^cm"^). We assume that the initial 



cloud evolution until a protostar is formed 
clouds have magnetic fields of Bq = 0-34 yuG, and angular velocities of Qq = 0-3.1 x 10~^"^s~^, 
as lis ted in Table [TJ These v alues are comparable with observed magnetic fields (ICrutcher 
19991 ) and angular velocities (iGoodman et al. Ill993l : ICasellill2002l ) of molecular cloud cores. 
Figure [2] shows the evolution of the magnetic field Be and angular velocity Qc at the center 
of clouds for Models 3, 8, 16, 22, and 23. The growth rates of the magnetic field for the non- 
ideal MHD Models 16, 22, and 23 (c^ ^ 0) are smaller than those for the ideal MHD Models 
3 and 8 (c^ ^ 0) for 10^^ cm~'^ ^ ""-c ^ 10^^ cm~^ (Fig. [2] upper panel). This is because the 
magnetic field is effectively dissipated by Ohmic dissipation for this phase. The magnetic 



fields at the star formation epoch (nc 



lO^^cm-^l 



B{ are listed in Table [TJ Models with 



= (ideal MHD models) have magnetic fields of ~ 100 kG (Model 3: 106 kG, Model 8: 
67.3 kG), while models with 1] (non-ideal MHD models) have magnetic fields of 0.1-1 kG 
(Model 16: 0.11 kG, Model 22: 0.68 kG, Model 23: 0.32 kG). The magnet ic field strengths of 
protos ta rs observed from Ze eman broadening measurements are ~ 1 kG (IJohns-KruU et al. 



1999al lbl.[ 200 ll : lBouvierll2006l ) . Even though protostars in non-ideal MHD models have masses 
of only ~ 10~^ Mq in our simulation, they have magnetic field strengths equivalent to those of 
observed protostars. The magnetic field is considered to be amplified by magneto-rotational 
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instability (MRI; iBalbus fc Hawleyi Il99ll ) and convection on the surface of a protostar in 
the later phase of star formation. However, our results indicate that protostars already 
have magnetic fields of ~kG for non-ideal MHD models at their formation epoch. On the 
other hand, when a protostar is formed, the magnetic field for ideal MHD models reaches 
B ~100 kG, which is about 100 times larger than the magnetic field of observed protostars. 

The rotation periods at the star formation epoch (n^ ~ 10^^ cm~^) P are listed in 
Table [TJ Figure [2] shows the evolutions of the angular velocity Qc (left axis) and rotation 
period P (right axis) at the centers of clouds. This figure shows that the rotation period 
reaches P = 1-100 days at the protostar formation epoch [ric = lO^^cm"^), even when the 
initial cloud has a small angular velocity (2 x 10~^''s~^). Many observations i ndicate that 



protostars have rotation periods of P = 0.6-20 days (e.g., iHerbst et al. II2006I ). Thus, the 
rotation periods derived from our calculation are roughly in agreement with those of observed 
protostars. The angular momentum is considered to be modified by interactions between 
the protostar and protoplanetary disk, jets, and outflow in the later phase of star formation. 
However, in our calculations, protostars have rotation rates of the same order of magnitude 
at their formation epoch, as is the case for the magnetic field. 

We show the evolution of the magnetic field in a collapsing cloud in the following sections 
in detail. Firstly, we show the cloud evolution for ideal MHD models [crj = 0) to simply 
investigate the evolution of the magnetic field. The cloud evolution for non-ideal MHD 
models (c^ ^ 0) is shown in subsequent sections. 



4.1. Evolution of the Magnetic Field for Ideal MHD Models 



4-1.1. Cloud Evolution with Different Parameters 
Figure [3] shows the evolution of the magnetic field normalized by the square root 



of the density p, 



1/2 



in the unit of i nitial sound speed c.^^n at the c enters of clouds for ideal 

1 /2 

0). As shown in iMachida et al. I (l2005al . I2006al ). the evolution of B^ pc 



MHD models (r/ 

depends on the mode of the cloud collapse, spherical, vertical, or disk-like collapse. When 
the cloud has a small amount of rotational or magnetic energy, the cloud collapses almost 
spherically (spherical collapse). On the other hand, the cloud collapses vertically along the 
rotation axis or magnetic field lines (vertical collapse) when the cloud has a large amount of 
rotational or magnetic energy. In either case, a thin disk finally forms around the center of 
a cloud, and the cloud continues to collapse keeping a disk-like structure (disk-like collapse). 
We briefly summarize the relatio n between the growth r ate of the central magnetic field and 
the cloud density (for details see lMachida et al. Il2005al . l2006al ). The magnetic field strength 
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increases as Be oc pl^^ in weakly magnetized and slowly rotating clouds (spherical collapse), 
because the cloud evolution is mainly controlled by thermal pressure and gravity. The 
magnetic field remains constant at oc p° when the cloud is magnetized strongly or rotating 
rapidly (vertical collapse), because the radial contraction of the cloud is suppressed by the 
strong magnetic tension or centrifugal force. After the disk formation, the magnetic field 

1/2 

evolves as 5c oc pc , and the Lorentz and centrifugal forces are balanced with the thermal 
pressure gradient and gravitational forces in the collapsing cloud (disk-like collapse jf]. For 
the isothermal phase (ric ^ 10^^ cm~^), -B//?^/^ at the center of a cloud converges to a certain 



value (the magnetic fiux-spin relation; iMachida et al. Il2005al ) when the cloud has a smaller 



rotational energy than the magnetic energy at the initial state. Thus, the initial strength 
of the magnetic field is not sensitive to the cloud evolution, especially for ric > 10^^ cm~^, 
if the field strength is sufficiently strong. Observations indicate th at the ma.gnetic energy 



is mu ch larger than the rotational energy in molecular cloud cores (jCrutcherl Il999l : ICaselli 
20021 ). Thus, convergence of B / p^^'^ is expected in the collapsing region of real molecular 



clouds. 

Figure [3] upper panel shows evolutions of the magnetic field in non-rotating clouds 
{u = 0) having different magnetic fields at the initial stage (Models 1, 2, 3, 4, 5, and 6). 
The magnetic field of Model 1 {Bq = 34 pG) is about 70 times stronger than that of Model 
6 {Bq = 0.5 pG) at the initial stage. The figure shows that the magnetic fields for all the 

1 /2 

models converge to certain values [BjpJ ~ 0.55 (Model 6) to 1.2 (Model 12) ] at the end 
of the isothermal phase {uc — 10^^ cm^^). Bc/pl^"^ begins to increase after the gas becomes 
adiabatic > lO^^cm^^), because the geometry of the collapse changes from disk-like to 
spherical for increasing thermal pressure. However, the clouds have almost the same magnetic 
field strengths for this epoch (nc > lO^^cm"^), because the growth rates of the magnetic 
field are almost the same for these models. The formed protostars have magnetic fields of 
Be = 88.9-123 kG at the protostar formation epoch [ric — lO^^cm"^), as listed in Table [TJ 
Thus, the formed protostars have magnetic field strengths different by a factor of 1.4, while 
clouds strength differ by a factor of 70 at the initial stage. As a result, protostars formed 
from initially non-rotating clouds have almost the same magnetic field if Ohmic dissipation 
is ignored. As shown in later sections, since both cloud rotation and Ohmic dissipation 
only decrease the magnetic field strength of a formed protostar, a protostar at its formation 
epoch cannot have a magnetic field exceeding B^ — 100 kG, derived from the non-rotating 
and ideal MHD models. 

The growth rate of the magnetic field for clouds with the same initial magnetic field 



In this paper, we use the terminology 'disk-hke collapse' when the disk is formed and the magnetic field 
increases as (x p'J'^ at the center of the cloud. 
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depends on the cloud rotation, because the geom etry of the collapse depends n ot only on 
the magnetic field but also on the cloud rotation (IMachida et al. Il2005al . l2006al ). Figure [3] 
lower panel shows the evolution of the magnetic field strength normalized by the square root 
of the density pd"^ for ideal MHD models (c^ = 0) with different initial angular velocities 
and the same initial magnetic field strengths (-Be = 1.6 /iG). The figure indicates that the 
growth rate of the magnetic field is smaller for an initially rapidly rotating cloud. Note that 
the calculation was stopped at n^. ~ 10^^ cm~^ in Model 11, with the largest angular velocity 
at the initial state, because fragmentation occurred. At the protostar formation epoch 
(nc = lO^^cm"'^), the cloud with the initially slowest angular velocity [lo = 0.001; Model 
7) had a magnetic field of = 95.9 kG, while the rapidly rotating cloud with uj = 0.1 
(Model 10) had Be = 43.7 kG (see Tabled]). This is because the growth rate of the magnetic 
field changes from B oc p^/^ to 5 oc p^/^ at an earlier evolutional stage in more rapidly 
rotating clouds. Therefore, the magnetic field at the protostar formation epoch depends 
on the rotation of the initial cloud. However, a slowly rotating cloud has a magnetic field 
strength only twice that of a rapidly rotating cloud (see Models 7-11 in Tabled]). 



4.1.2. Cloud Evolution of Typical Model 



Figure m shows the cloud evolution from the initial state [ric = 1.4 x lO^cm"'^; Fig. Hja) 
to the protostar formation epoch [ric = ^.^ x 10^^ cm~^; Fig.H]) for Model 12 with parameters 
a = 0.01, UJ = 0.01, and = 0. In this model, the magnetic field is well coupled to the gas 
from the initial to the final stage because c,, = 0. The initial cloud is weakly magnetized 
(-Bo = 1.4 pG) and rotating slowly {Qq = 2 x 1 0~^^ s~^). This cloud is classified as a 



magnetic-force dominant model by the criterion in iMachida et al. I (j2005al . l2006al ). because 
the ratio of the initial angular velocity Qq to the magnetic field Bq is i^oZ-^o = 1.4 x 10^^ < 
5.3 X 10~^ = {Q/B)cri- Thus, the magnetic field affects the cloud e volutio n, while the effect 
of cloud rotation is small (for details, see IMachida et al. I l2005a] . l2006al ). Figure [5] shows 
the magnetic field distribution at the same epoch for the panels in Figure HI As shown in 
Figure Hja and b, the central region gradually becomes oblate as the cloud collapses and a thin 
disk is formed in the isothermal phase {uc ^ 10^^ cm~^). The magnetic field lines gradually 
converge towards the center (F ig. [5]6 and c). This configuration of the magnetic field lines 
is similar to that observed by iGirart et al. I (120061 ). The thin solid lines in Figure [6] left 
panels show the evolution of the magnetic field normalized by the square root of the density 
(-Bc/p^^; Fig.Ejfl) in the unit of initial sound speed Cg^o, plasma beta {f3p = Sirclp/ B'^; Fig.[H]&) 
within the region p > O.lpc, and specific angular momentum normalized by the mass [J/M"^; 
Fig. [6]c) within the region p > O.lpc in the unit of AnG/csfi, against the central density. 
Figure [6]a and b show that both Bc/pl^"^ and J/M'^ continue to increase in the isothermal 
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phase {ric < 10 cm~ ). Thus, it is considered that the disk is formed by the Lorentz and 
centrifugal forces, as shown in Figure [5]c. The plasma beta has a minimum of /9p ~ 8 at the 
end of the isothermal phase {n^ — 10^^ cm~^ — 10^^ cm~^; Fig.lHlmiddle panel). Even after the 
gas around the center of the cloud becomes adiabatic {ric > lO^^cm"^), the magnetic field 
continues to be amplified (Fig. [6]a). On the other hand, the plasma beta remains (3p ~ 10-20 
for He ^ lO^^cm"^, after it increases slightly owing to an increase of the thermal energy for 



10 



12 



cm 



-3 < 



< 10 



16 



cm 



When the gas reaches ric — 10 cm" , the central region becomes optically thick and 
the equation of state b ecomes hard, as derived from the one-dimensional radiative hydro- 
dynamical calculation (jMasunaga &: Inutsuka I l2000l ) . After the equation of state becomes 



hard, shock occurs and the first core is formed. The first core is formed at ~ 10^'^ cm~'^ 
in this model. The first core is shown by the white-dotted line in Figure Hlc and d. (We plot 
the position where the shock occurred on the z = plane in these panels.) The first core has 
a mass of ~ 0.012 and radius ~ 7AU. The arrows in Figure indicate that the radial 
component of the velocity fr is dominant outside the first core, while the azimuthal com- 
ponent of the velocity is comparable to the radial component inside the first core on the 
z = plane. This is because the large thermal pressure suppresses the cloud collapse inside 
the first core, and thus the rotation timescale becomes smaller than the collapse timescale. 
Figure lower panel shows that the outflow is driven from the first coreE Outflow appears 
when the centr al density reaches rir. — 10^^ cm~^. Outflow driven fr om the first core has als o 



been shown bv iToinisakal (120021) iMatsumoto fc Tomisakal (120041 ). iMachida et al. I (l2005bf ). 
Banerjee &: Pudrita (120061 ) , and iFromang et al. I (120061 ) . As shown in Figure EH, the mag- 
netic field lines begin to twist due to the rotation of the first core. This outflow reaches 
~ 10 AU and has a maximum speed of ~ 3kms~^ at the end of the calculation. 

The specific angular momenturn norma lized by the mass J/M^ begins to decrease in the 
adiabatic phase. IMatsumoto et al. I (119971 ) and ISaigo fc Tomisakal (120061 ) show that J/M^ 
is constant after it reaches a peak in unmagnetized clouds. Thus, the decrease of J/M^ is 
caused by magnetic effects. Since no outflow appears before the first core formation, th e 
decrease of J/M^ for this phase is caused by magnetic braking (IBasu fc Mouschoviaslll994l ). 



After the central density reaches Uc — 10 cm~ , J/M decreases more rapidly. This shows 
that the angular mor nentum is remov ed n ot only by magnetic braking, bu t also by outflow in 
the adiabatic phase. iTomisakal (119981 ) and IMatsumoto fc Tomisakal (120041 ) show that outflow 
is an important mechanism for angular momentum transfer. 



The cloud collapses again inside the first core for 10 cm 



,-3 < 



rir 



< 10 cm because 



We call the flow driven from the first core 'outflow' and the flow driven from the second core a 'jet'. 
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of dissociation of molecular hydrogen (i.e., the second collapse). Then, the equation of state 
becomes hard again and the second core (i.e., a protostar) is formed at n^. ^ 10^^ cm~^ 
(IMasunaga fc Inutsuka II2OOOI ). Figure Hie shows the structure near the protostar (the proto- 
star or shocked region is represented by the black-dotted line). The gas around the protostar 
has a density range of 10^^ cm~^ ^ ''^c ^ 10^° cm~^ and collapses rapidly, because the equa- 
tion of state is soft in this region [see Eq. ([5])]. Thus, the gas accretes onto the protostar 
with a high speed of ~ lOkms"^. Since the rotation velocity is comparable to the accretion 
velocity, the magnetic field lines are strongly twisted (Fig. [5](i). Figure ll]f shows the structure 
around the protostar 215 hours after the protostar formation epoch. The protostar has a 
mass of 2.1 x 10"'^ and a radius of 1.1 i?© at this epoch. The disk surrounding the pro- 
tostar extends up to 3.9 i?© with 1.2 x 1O~^M0. A strong jet is driven from the protostar, 
as shown i n the u pper panel of Figure BIf. Th e jet driven from the protosta r is also shown by 
Tomisakal J2002h . iBaneriee fc Pudr"it3T2006h . and lMachida et al. I J2006bh . The jet reaches 



lli?0 at the end of the calculation. In this model, since the effect of Ohmic dissipation is 



ignored 



0), the magnetic field couples with the gas at any time. Thus, the magnetic 



field around the protostar is strongly twisted at the protostar form ation epoch, as s hown 
in Figure Elf- The c l oud e volution of Model 12 is similar to that of iTomisakal (120021 ) and 



Banerjee fc Pudrita (l2006l ). in which the evolution of a magnetized cloud in an ideal MHD 



regime is studied. 



4.2. Evolution of Magnetic Field in Non-Ideal MHD Models 

Figures [71 and [H show the cloud evolution for Model 16. Model 16 has the same magnetic 
field strength (a = 0.01, Bq = 1.4 ;uG) and angular velocity {u = 0.01, Qq = 1.4 x 10~^^s"^) 
as Model 12 (Figs. HI and ED at the initial state. However, Model 16 has = 1, which means 
that the magnetic field is dissipated by Ohmic dissipation in a collapsing cloud (non-ideal 
MHD model), while the effect of Ohmic dissipation is ignored in Model 12 [crj = 0; ideal 
MHD model). In Model 16, the magnetic Reynolds number is Re^ < 1 for 2 x lO^^cm"^ < 
nc < 6 X 10^^ cm~^ (Fig- [ID- Thus, the cloud evolutions of Models 12 and 16 are almost 
the same for Uc <^ 10^^ cm"'^, because Ohmic dissipation is not effective for this early phase 
even in Model 16. The magnetic field normalized by the square root of the density B/p^^'^, 
plasma beta /3p, and specific angular momentum normalized by mass J/M"^ in Model 16 
(thick solid lines in Fig. [Sla-c) are identical to those in Model 12 (thin solid lines in Fig. [Hla- 
c) for Tic ^ 10^^ cm~^. Thus, Figures [3 and [HI show the cloud evolution for Model 16 only 
after the first core is formed. 

The first core is formed at ric — 10^^ cm~^ for Model 16, as for Model 12. Figures [71a and 
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HJc show that the first cores in Models 12 and 16 have almost the same shape and size. The 
first core has a mass of 0.012 Mq and a radius of 7AU in Model 16. As shown in Figure [8]a, 
the poloidal component of the magnetic field is dominant at the first core formation epoch 
[ric — lO^^cm"^). After the formation of the first core, the magnetic field lines begin to be 
twisted (Fig.P). 

The outflow from the first core is shown in the lower panel of Figure [7]c (red solid lines) . 
Outflow appears 22 yr after the first core formation epoch in Model 16 {uc — 10^''cm~^), 
while outflow appears 2.6 yr after the first core formation epoch {ric ~ 10^^ cm~^) in Model 
12. In Model 16, the outflow has a maximum speed of ~ 2kms~^ and reaches ~ 1 AU at the 
end of the calculation. Figures [7|c and show that the outflow in Model 16 is smaller and 
slower than that in Model 12. The outflow is mainly driven by the Lorentz force outside the 
first core. The magnetic field is weak and barely twisted around the first core in Model 16, 
because the magnetic field in this region is dissipated by Ohmic dissipation. Thus, outflow 
in non-ideal MHD models (Model 16) is weaker than that in ideal MHD models (model 12). 
We will further discuss the driving mechanism of outflow and jets in a companion paper. 

Figure E]a shows that the magnetic field in Model 16 (thick solid line) is largely removed 
from the central region for lO^^cm"'^ ^ '^c ^ lO^^cm"^. The plasma beta in Model 16 
increases from /5p ^ 10 to ~ 2 x 10^ for the same epoch (Fig. [6]&). In Figure [8]c, the magnetic 
field lines are strongly twisted near the center of the cloud, while they are barely twisted 
away from the center of the cloud. This implies that the neutral gas is well coupled with 
the magnetic field, and the magnetic field lines rotate with the neutral gas near the center 
because Ohmic dissipation is ineffective. When the gas density reaches ric — 10^^ cm~^, 
the temperature reaches T ~ 2000K. Thus, the thermal ionization of alkali metals reduces 
the resistivity and Ohmic dissipation becomes ineffective. As a result, the magnetic field is 
coupled with the gas and amplified again for Uc ^ 10^^ cm~^. On the other hand, the neutral 
gas is decoupled from the magnetic field in the region away from the center and the magnetic 
field lines slip through the neutral gas and rotate freely because this region has a density range 
of 10^^ cm~^ ^ '^c ^ 10^^ cm~^ (Fig- Wp lower panel) and magnetic dissipation is effective. 
Thus, the magnetic field lines that are strongly twisted in the earlier phase {uc -C 10^^ cm~^) 
are de-coUimated and relaxed by the magnetic tension force. Note that the magnetic field 
lines do not disappear, even with Ohmic dissipation. In our settings, it is assumed the 
magnetic field has a vertical component in the initial state. Thus, even when both the 
azimuthal and radial B,- components of the magnetic field are completely dissipated 
by Ohmic dissipation, the vertical component of the magnetic field remains because of the 
condition V ■ -B = 0. 

As shown in the upper panel of Figure [7H, the gas inside the first core rotates rapidly. 
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The specific angular momentum in Model 16 is larger than that in Model 12 (Fig. [6]c) 
for rtc > lO^^cm"^, while the magnetic field in Model 16 is smaller than that in Model 12 
(Fig.[H]a)- This means that magnetic braking is less effective in Model 16. Therefore, a strong 
centrifugal force creates a thin disk at the center of the cloud, as shown by the projected 
density contour on the sidewall in F igure and e. In additi on, the cloud collapses slowly 



for a rapid rotation, as discussed in ISaigo fc Tomisakal (120061 ). It takes 40.7 yr to form the 



second core {uc — 10^^ cm~^) after the first core formation epoch (rzc — 10^^ cm~^) in Model 
16, while it takes 25.5 yr in Model 12. The magnetic field lines become twisted again inside 
the first core (Fig. [8](i and e). Then, a different flow from that driven from the first core 
appears, as shown in the lower panels in Figures and e. Inside the red line in Figure [7](i 
lower panel, slow flow (outflow) driven from the first core is seen in the region \z\ > 0.4, 
while fast flow (a jet) driven from a rotational supported core inside the first core is seen in 
the region \z\ < 0.4. 

To investigate the properties of the outflow and jets, we show in Figure [9] the z- 
component of the velocity Vz (upper side) and the plasma beta Pp (lower side) at the same 
epoch and with the same scale as Figure [7]d. In this figure, we can see two sets of peaks 
in the 2;-component of the velocity (outer weak and inner strong peaks; Fig. [H] upper side). 
The outer weak peaks with Vz — 5kms~^ are located at (x, z) = (±0.3, 0.3) and originate 
from the first core. The inner strong peaks with Vz — 15kms~^ are located at (x, z) = 
(±0.05, 0.3) and originate from the rotation supported quasi-static core inside the first core. 
In the upper panel, we label the flow driven from the first core as 'outflow', and the flow 
driven from the quasi-static core as 'jet'. Figure [3 lower panel shows that the plasma beta 
is /9p < 1 around the velocity peak, which indicates that flows with t>^ > or t>2 < for 
z > OT z < 0, respectively, are driven from the region with strong magnetic field. These 
regions are composed of low-density gas {uc ^ lO^^cm"^), as shown in Figure [7](i, and thus 
Ohmic dissipation is ineffective (Fig. [T]) and the magnetic field can be amplified as the cloud 
collapses. On the other hand, the central object and disk located at z < |0.1| have a large 
plasma beta [jSp = 10^ — 10^) because Ohmic dissipation is effective in this region. The jet is 
driven from the second core (protostar) after the protostar formation epoch [ric = 10^^ cm^^) 
in Model 12, while the jet is driven from the quasi-static core before the protostar formation 
epoch in Model 16. This is because the rotation supported core is formed before the central 
density reaches Uc = 10^^ cm~^ in Model 16, due to lower effectiveness of magnetic braking. 
Figure [7]e shows the cloud structure when the density reaches nc = 8.1 x lO^^cm"^. In this 
figure, a strong jet is driven from the center. 

When the gas density reaches ric ~ 1.3 x 10^^ cm~^ in Model 16, a shock occurs near the 
center of the cloud and the protostar (second core) is formed. The protostar in Model 16 is 
formed at a lower density than that in Model 12. The protostar has a mass of 9.1 x 10~^ Mq 
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and a radius of 3Rq at its formation epoch. The radius of the protostar in Model 16 is 
three times larger than that in Model 12. The central region inside the first core rotates 
rapidly in Model 16, thus, a considerably flattened disk (Fig. [3f lower panel) is formed near 
the center due to the strong centrifugal force. Figures U]f and ^ show that the protostar in 
Model 16 rotates rapidly and has an oblate structure. Figure [8]f shows that the magnetic 
field lines are weakly twisted around the protostar, where Ohmic dissipation is effective. 
This is because the rapid rotation makes the first core (or high-density region) increasingly 
oblate as the cloud evolves, and the regions above and below the first core and protostar have 
low-density gases [ric < lO^^cm"^), as shown in Figure [7]f lower panel. Thus, the magnetic 
field is decoupled from the neutral gas and again de-coUimated in this region. 

As a result, Ohmic dissipation affects not only the distribution and strength of the 
magnetic field but also the cloud collapse, rotation, outflow, and jet. 



4.3. Dependence on Cr, 

Ohmic dissipation greatly influences the cloud evolution, as described in the previous 
se ction. We estirnate th e resistivity as a function of the density and temperature according 
to iNakano et al. I (120021 ) . shown by the thick solid line in Figured] However, we note that 
it is difficult to determine the ionization structure accurately owing to the uncertainties, 
particularly in the properties of dust grains that provide the most important sites for recom- 
bination in collapsing clouds. Thus, we parameterize the model of the magnetic dissipation 
Crj and investigate the cloud evolution, as shown in Figure [61 Figure [6] left panels show the 
evolution of models with the same magnetic field and angular velocity {a = 0.01, u = 0.01) 
but different c^. In Figure Eh, the larger efficiently removes the magnetic field from the 
cloud. The plasma beta is maintained at (3p ~ 10-20 for Hc > 10^^ cm~^ in Model 12 (the 
ideal MHD model; = 0), while it reaches /3p ~ 10^ for 10^^ ^ "-c ^ 10^^ cm"^ in Model 16 
[Cfj = 1) in Figure [6f). Thus, the central region in Model 12 has magnetic energy ~ 10^ times 
larger than that in Model 16. This difference in the magnetic energy between Models 12 and 
16 is caused by Ohmic dissipation for lO^^cm"^ ^ ^ lO^^cm"'^. The magnetic energy in 
the central region depends on Cr^, as shown in Figure [6]6. In contrast, the angular momentum 
is large and the cloud rotates rapidly in models with larg shown in Figure [6]c. This is 

because magnetic braking is less effective in models with larger for less magnetic energy. 
As a result, protostars with larger have larger angular momentum at the formation epoch. 
In Model 17 {crj = 10), the cloud rotates rapidly and fragmentation occurs at ric — 10^^ cm~^, 
because the magnetic energy is small and magnetic braking is barely effective. 



Figure E] right panels (Fig. (HH-/) show the evolution of models having the same magnetic 



- 16 - 



field [a = 0.01) and rotation rate [uj = 10~^) with different c^. In the initial state, the 
rotation rates of the models shown in Figure [6] right panels are 100 times smaller than those 
in the left panels {uj = 10~^), while the magnetic field strengths are the same. The evolution 
of the magnetic field and plasma beta is almost the same for ric ^ 10^^ cm~'^ in both FigureEh 
and d. After the magnetic dissipation, the evolutions of the models shown in the left and 
right panels are different. The value Be/ pd"^ increases with gas density for ric ^ 10^^ cm~^ in 
Models 18-23 (Fig. [6]d), while it is almost constant in Models Vl-Yl (Fig. [6]a). The magnetic 
field is more greatly amplified for a spherical collapse [B oc p^^^) than for a disk-like collapse 
[B oc p^^^), as explained in §4.1. In ideal MHD models, even if the cloud has no angular 
momentum, the geometry of the collapse changes from spherical to disk-like collapse, because 
the magnetic field is amplified and a disk-like structure is formed at the center of the cloud. 
After the disk-like structure is form ed (i.e., during the disk -like collapse), the growth rate 
of the magnetic field becomes small (IMachida et al. Il2005al ). However, if the magnetic field 
is removed from the central region sufficiently, the geometry of the collapse again changes 
from disk-like to spherical collapse and the magnetic field is amplified as S oc (spherical 
collapse). In Models 18-21, after the dissipation of the magnetic field (n^ > lO^^cm"^), the 
clouds collapse spherically because of the extremely slow rotation, while in Models 12-17 the 
clouds collapse with a disk-like shape due to the rapid rotation. At the protostar formation 
epoch (nc ~ 10^^ cm~^), the magnetic field is amplified to i?c = 0.11 kG in Model 16 (rapidly 
rotating cloud), while it is amplified to B^ = 0.68 kG in Model 22 (slowly rotating cloud). 
Thus, the protostar has a stronger magnetic field for a more slowly rotating cloud, even in 
non-ideal MHD models. After the magnetic field is significantly removed from the center, the 
angular velocity is also amplified because the cloud collapses spherically. Thus, the specific 
angular momentum normalized by the mass J/M^ also increases in Models 18-23 (Fig. [6]f), 
while J/M^ is constant in the more rapidly rotating clouds of Models 12-17 for disk-like 
collapse (Fig. [6]c). 



4.4. Cloud Evolution with Different Angular Momentum 



Figure [10] shows the evolution of the magnetic field Be/ pc (upper panel) and angular 
velocity f2/(47rGpc)^^^ (lower panel) for Models 24-28. These models have the same mag- 
netic field (a = 0.01) and (c^ = 1), but different angular velocities u. In this figure, 
the circles indicate the fragmentation epoch. Fra gmentation occu r s in th e rapidly rotating 



clouds of Models 27 {u = 0.03) and 28 {uj = 0.1). iMachida et al. I fl2005bf ) shows that frag- 



mentation occurs when the central angular velocity Qc normalized by the free-fall timescale 
l/(47rGpc)^/^ reaches f2c/(47rGp)^/^ ~ 0.2 for the isothermal phase ric ^ lO^^cm"^. As 
shown in Figure [10] lower panel, this condition is realized in Models 27 and 28. After frag- 
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mentation occurs, we stopped the calculation because our numerical code (nested grid) is 
not suitable for calculating the evolution of each fragment located outside the center. The 
protostar forms without fragmentation in Models 24, 25, and 26, because the fragmentation 
condition is not realized in these clouds, as shown in Figure dU] lower panel. 

Figure [TO] upper panel shows the evolution of the magnetic field B^/ pc ■ In this panel, 
the clouds without fragmentation have almost the same magnetic field strength for < 
lO^^cm"'^. For ric ^ lO^^cm"^, the magnetic field continues to increase in the non-rotating 
cloud (Model 24), while they saturate at certain values in the rotating clouds (Models 25 
and 26). A close analysis of the evolutions for ric > 10^^ cm~^ shows that the growth 
rate of the magnetic field becomes small when the normalized angular velocity approaches 
Vic/ (4vrGpc)^''^ — 0.2. This is due to the fact that the collapse of the geometry changes from 
spherical to disk-like when the rotational ene rgy becomes comparable to the gravitational 
energy [for details, see lMachida et al. Il2005a| . The protostars have magnetic fields of Be = 
0.51 kG (Model 24; uj = 0), 0.48 kG (Model 25; cu = 0.001), and 0.11 kG (Model 26; u = 0.01). 
In this way, the magnetic field strengths of formed protostars are related to the initial cloud 
rotation. In summary, a protostar has a stronger magnetic field in a more slowly rotating 
cloud. 



Observation s indicate that molecular clouds have larg e amounts of magnetic energy 
(ICrutcherl Il999l ) and small amounts of rotational energy (iGoodman et al. I Il993l : ICaselli 
2OO2I ). These clouds are classifie d as magnetic- force dominant clouds {flo/Bo < Vtcn/ B^n = 
0.39 G^/^ cJq.) in the criterion of iMachida et al. I (2005aj, |2006aj). As shown in §4.1, the mag- 
netic field converges to a certain value for the isothermal phase in magnetic-force dominant 
clouds. After the magnetic dissipation, the growth rate of the magnetic field depends on the 
rotation rate at the center of the cloud. Thus, the magnetic field of a protostar is determined 
almost completely by the rotation energy of the initial cloud when the degree of resistivity 



is fixed. 



5. Discussion 

5.1. Magnetic Flux Problem 

The magnetic flux problem is important for the star formation process. Magnetic energy 
in an interstellar cloud is usually regarded as being comparable to the gravitational energy of 
the cloud, while the magnetic energy in a star is much smaller than the gravitational energy. 
The molecular cloud can collapse when the following condition is realized: 



(12) 
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where $ is the magnetic flux and is defined as 



B-dS. 



and •I'cri is the critical value of the collapse and is defined as 



(13) 



(14) 



wher e / is a dimensionless constant (/ ~ 8; iMouschoyias fc Spitzerlll976l : iTomisaka et al. 
1998 ) and M is the mass o f the cloud or cloud cor^. The "observed magnetic fiuxes of 



molecular cloud cores (e.g., ICrutcherl Il999l ) are close to their critical value ($/$„! ~ 1), 
while those of protostar s are much smaller than their critical value ($/$cri ~ 10"^-10^^; 



Johns-KruU et al. 11200 ll ). Thus, the magnetic fiux should be largely removed from the cloud 
core during the star formation process. This is called "the magnetic fiux problem" in star 
formation. 

The fiux to mass ratio $/M aro und the center of a collapsing cloud de creases by two 
mechanisms. One is v ertical collapse (INakano Ill983l : iMachida et al. Il2005al ) and the other 
is Ohmic dissipation (INakano et al. 1120021 ). When a gas infalls vertically along the mag- 
netic field line (vertical collapse), the mass around the central region increases, keeping the 
magnetic fiux constant. Thus, the ratio of the magnetic fiux to the critical mass, $/$cri 
(cx $/M), decreases, beca use the denomin ator of $/$cri increases, while the numerator is 
nearly constant. However, INakano I (119831 ) showed that the vertical collapse mechanism is 
not sufficient for solving the magnetic fiux problem, because a large amount of gas collects 
from a wide range (> 200 kpc) to d ecrease this ratio from $/$cri — 1 to ~ 10^^-10"'^. On the 
other hand, INakano et al. I (120021 ) showed analytically that the Ohmic dissipation process 
in a collapsing cloud is an effective mechanism for solving the magnetic flux problem. We 
numerically conflrmed that Ohmic dissipation is more important than vertical collapse for 
the removal (or decrease) of magnetic flux in a collapsing cloud. 

Figure [TTb shows the distribution of the magnetic flux normalized by the critical mass 
$/$cri as a function of the cumulative mass from the center of the cloud at the initial state 
(broken line), and at the end of the calculation for Models 29 (solid line) and 30 (dotted 
line). In this graph, the mass M (x-axis) is integrated for every iso-density contour pa from 



^ In general, <& and <i>cri are used to measure the degree of the magnetic field strength of an isolated cloud 
core. In this paper, however, we use these values to measure the degree of the magnetic field strength inside 
a collapsing cloud core. We define <& and <I>cri simply as functions of the magnetic field and mass in the local 
region inside the cloud core. 
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the center that has maximum density, as 

M(p>pa) = / pdV, (15) 

and the magnetic flux inside the corresponding region is integrated as 

^{p >P.)= [ B dS, (16) 

where dS is defined at the z = plane and parallel to the ^-axis. The graph indicates 
that Models 29 and 30 have magnetically supercritical clouds {^/^cii < 1) as a whole 
(M > 0.3 Mq), while the central regions (M < 0.3 M©) of these clouds are magnetically 
subcritical ($/$cri < 1) at the initial state. The ratio ^/^rr\ increases with decre asing 



mass. Since we adopt the Bonnor-Ebert density profile (iBonnorl Il956l : lEbertl Il955l ). the 
initial density distribution near the center is almost constant and thus the mass of the cloud 
(i.e., $cri) is proportional to R^, M oc $cri oc R^. The magnetic fiux is proportional to 
$ oc i?^, because a uniform magnetic field is assumed at the initial state. Thus, the central 
region (comprising a small fraction of the total mass) has large $/$cri; since the ratio of the 
magnetic fiux to the critical mass is proportional to $/$cri oc R^^ oc M^^/^ near the center. 
However, these clouds can collapse promptly, because a large amount of the mass (~ 99% of 
the total mass) is magnetically supercritical. 

In our calculations, protostars have a mass of ~1O^'^M0 at their formation epoch. 
Figure [Tib shows that the ratio of the magnetic fiux to the critical mass within a mass 
of M = 10"^ M0 (i.e., at M = 10"^ on the x-axis) decreases from $/$cri = 6 to 
$/$cri — 0.2 in Model 30 (ideal MHD model). In Model 30, the ratio $/$cri decreases due 
to vertical collapse, because Ohmic dissipation is ignored since c,, = 0. This indic ates that 



the m agnetic fiux is not removed sufficiently by the vertical collapse, as shown in iNakano 



(119831 ). On the other hand, in Model 29 (non-ideal MHD model), the ratio $/$cri within 
M = 10"^ M0 decreases from $/$cri = 6 to <I>/<I>cri — 10"^ (Fig. [TTja). Thus, three or 
four orders of magnitude of the initial magnetic fiux are removed from the core by Ohmic 
dissipation in Model 29. The solid line in Figure [TTb shows that the ratio $/$cri decreases 
from $/$cri ^ 0.1 to <l>/$cri ^ 10~^ steeply between lO-^M© < M < 5 x lO'^M©. This 
drop indicates that the gas within M < 10~^ Mq from the center of the cloud experiences 
Ohmic dissipation. The remainder of the cloud is expected to experience Ohmic dissipation 
as the collapse continues. 

Figure [TTb shows the evolution of $/$cri for Models 29-32 against the central density. 
In this graph, the magnetic fiux $ and the mass M within the gas having density p > O.lpc 
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are integrated at each timestep as 

M(p>0.1pc)= / pdV, (17) 

<l>(p>0.1pc) = / B-dS. (18) 

J/3>0.1pc,Z = 

^/'I'cri in Models 30 and 32 (ideal MHD models) decreases slightly from the initial state 
($/$cri — 0.8) to the protostar formation epoch ($/$cri — 0.2). In contrast, for Models 29 
and 31 (non-ideal MHD models), $/$cri decreases significantly for 2 x 10^^ cm~^ ^ '^c ^ 
6 X 10^^ cm~^ by Ohmic dissipation, and reaches $/$cri — 10"'^ at the protostar formation 
epoch. 

In the non-ideal MHD models, the formed protostars have smaller $/$cri ratios for 
the more rapidly rotating clouds. After the magnetic dissipation (ric > lO^^cm"^), the 
evolution of $/$cri is different in the non-ideal Models 29 and 31 (Fig. [TTb). For this phase 
{ric > lO^^cm"'^), the ratio is almost constant ($/$cri — 5 x 10~^) in Model 31 (with a 
rapidly rotating cloud), while it continues to increase in Model 29 (with a slowly rotating 
cloud). At the protostar formation epoch, the protostars have $/$cri = 2 x 10"'^ for Model 
29 and $/$cri = 4 x 10~^ for Model 31. This difference is caused by the rotation rate of 
the initial cloud. In Model 31, the cloud collapses vertically (i.e., vertical collapse) and the 
growth rate of the magnetic field is low because of the rapid rotation for Uc ^ 10^^ cm~^. On 
the other hand, in Model 29, the cloud collapses spherically (i.e., spherical collapse) and the 
growth rate of the magnetic field is high for ric ^ 10^^ cm~^, because both the Lorentz and 
centrifugal forces are weak. After the magnetic dissipation {ric > 10^^), since the Lorentz 
force is weak, the geometry of the collapse is determined only by the cloud rotation. An 
initially slowly rotating cloud collapses spherically to form a protostar with a large magnetic 
flux, while a rapidly rotating cloud collapses along the rotation axis to form a protostar with 
a small magnetic flux. 

In summary, protostars have $/$cri — 10"'^ to 10~^ at their formation epoch, which 
is comparable to observations. However, since we stopped our calculation at Mprotostar ~ 
1O~^M0, further calculations are needed to predict the magnetic flux for protostars with 

Mprotostar ~ 1 Mq. 

5.2. Angular Momentum Problem 

The angular momentum problem is also important for the star formation process. 
Molecular cloud cores have specific angular momenta of the order of jcioud ~ lO^^cm^ s~^ 
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( IGoodman et al. I Il993l : ICasellil |2002| ). while protostars have specific angular momenta of 
the order of jprotostar ~ lO^^cm^s"^. Thus, the specific angular momentum of a protostar 
is ~ 10~^ times smaller than that of a molecular cloud core. For this reason, the angular 
momentum of the initial cloud must be removed from the core as well as the magnetic flux. 
Tomisaka J2OO0I . [ 2OO2I ) calculated the cloud evolution from a molecular cloud to protostar 
formation paying attention to the evolution of the angular momentum in two-dimensional 
nested grid simulations. The study showed that the angular momentum of the initial cloud is 
sufficiently removed by magnetic braking and outflow in the collapsing cloud, and the formed 
protostar has a specific angular momentum comparable to observed values. However, the 
ideal MHD approximation was used and thus the magnetic effect (i.e., magnetic braking 
and outfiow) may have been overestimated. In this section, we discuss the evolution of the 
angular momentum in the non-ideal MHD regime. 

Figure [TTb shows the distributions of the specific angular momentum against the cumu- 
lative mass at the initial state (thick broken line), and the end of the calculation for Models 
29 (solid line), 30 (thin broken line), and 33 (dash-dotted lines). In this figure, the specific 
angular momentum is integrated up to every iso-density surface pa from the center as 



j{p > Pa) 



M(p > Pa 



/ pzu -v^dV, 

J P>P3. 



(19) 



where M(p > pa) is the mass within p > Pa [see Eq. (IT7|) ] and w is the radius in cylindrical 
coordinates. As shown by the thick dashed line, the specific angular momentum at the initial 
state is proportional to jcioud oc near the center of the cloud, starting with rigid rotation 
and the Bonnor-Ebert density profile with an almost constant density at the center. 

In Figure [TTb . the specific angular momentum at the end of the calculation (jprotostar) 
for Model 33 is shown as a dash-dotted line. In Model 33, the specific angular momentum 
is not removed by magnetic effects (magnetic braking, outflow and jet), because the cloud 
is unmagnetized. Thus, the distribution of the speciflc angular momentum at the end of the 
calculation is almost the same as that of the initial state for M > 10^^ M©. This indicates 
that the angular momentum is almost conserved outside the protostar (M > IO^^Mq). 
For M < 1O~^M0, the speciflc angular momentum at the end of the calculation is slightly 
smaller than that at the initial state, caused by angular momentum transfer by the non- 
axisymmetric pattern. In this model, after a thin disk is formed, the disk is transformed 
into a bar, because bar mode instability (e.g., iDurisen et al. Ill986l ) is caused by rapid ro- 



tation. The speciflc angular momentum begins to decrease around the center of the cloud 
as the non-axisymmetric pattern grows. Fragm entation in this unmagnet ized model is ob- 
served for further calcul ations, as also shown by lBanerjee fc Pudritzl (120061 ) . As discussed in 
Machida et al. I (j2005bl ). fragmentation plays an important role in the angular momentum 
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evolution, because the angular momentum is redistributed into the orbital and spin angular 
momentum of each fragment. We will discuss the fragmentation and redistribution of the 
angular momentum in subsequent papers. 

The specific angular momenta of magnetized clouds at the end of the calculation are 
plotted by the solid (non-ideal MHD Model 29) and thin broken lines (ideal MHD Model 
30) in Figure [Tib. The specific angular momenta of both Models 29 and 30 at the end of the 
calculation are smaller than those at the initial state. In these models, the specific angular 
momentum decreases by vertical collapse, magnetic braking, outflow, and jets. When the 
initial cloud is magnetized strongly (or rotates rapidly), the cloud collapses vertically along 
the magnetic field line (i.e., vertical collapse). Gas with small specific angular momentum 
preferentially falls into the central region, and thus the central region gains specific angular 
momentum in the vertical collapse regime. Note that vertical collapse does not transfer 
angular momentum, because this mechanism changes only the location of each fluid element 
and thus the gas keeps its specific angular momentum. For higher densities {ric > 10^^ cm~^), 
the magnetic field lines are strongly twisted, and outflow and jets appear. T he angular 



momentum is effectively transferred by the outflow and jets, as already shown by iTomisaka 



( I2OOOI . I2OO2I ). Magnetic braking is also important for the angular momentum transfer. In 
principle, however, a reduction of the angular momentum due to magnetic braking cannot be 
distinguished from that due to outflow and jets, because outflow and jets are consequences 
of the magnetic braking process. 

Figure [TTl fi shows the evolution of the specific angular momentum j in the inner region 
for Models 29, 30, and 33 as a function of the central density. The angular momentum is 
integrated within p > O.lpc at each timestep and averaged as 



j(p > O.lpc) = , [ pw-v^dV, (20) 



M{p > O.lpc 



where M{p > O.lpc) is the mass within p > O.lpc- In this graph, the specific angular momenta 
in the magnetized clouds (Models 29 and 30) are smaller than that in the unmagnetized 
cloud (Model 33) for ric > 10^ cm^'^. The difference of the specific angular momentum 
for the isothermal phase {n^ < 10^^ cm~^) is due to magnetic braking and vertical collapse, 
because no outflow appears in this phase. The specific angular momentum in Model 30 (ideal 
MHD model) is smaller than that in Model 29 (non-ideal MHD model) for Uc > 10^^ cm~^, 
because magnetic braking is less effective due to magnetic dissipation for 2 x 10^^ cm~^ < 
^ 6 X 10^^ cm"^. 

At the protostar formation epoch [ric = lO^^cm"^), protostars have specific angular 
momenta of jprotostar = 8 x 10^^ cm^ s~^ for Model 31, 1.5 x lO^^cm^s"^ for Model 29, and 
2 X 10^^ cm^ s~^ for Model 30. Observations indicate that classical T Tauri stars have specific 
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angu lar momenta of the order of jobs ~ lO^^cm^s ^ (iBouvier et al. 1 119931 : iHerbst et al 



20061 ) ■ At the end of our calculation, the formed protostars have masses of only ~ 10~^ Mq. 
Since the spatial distribution of the specific angular momentum follows j oc cu^ for rigid 
rotation, the region inside M < 10~^ Mq of T Tauri stars has an average specific angular 
momentum of ~10^^cm^ s~^. Thus, the specific angular momenta in non-ideal MHD models 
such as Model 29 are about one order of magnitude larger than those of observed protostars. 
In this paper, we discuss the angular momentum evolution only in the early star formation 
phase (Mprotostar < 10^^ — 10"^ Mq). In the later accretion phase, however, the angular 
momentum is transferred by star-disk interaction, stellar wind, and jets. To fully solve the 
angular momentum problem, the evolution in the later phase must also be calculated. 

This paper focuses on the formation of sin gle stars, and hence we ha v e investigated the 

evolut ion of slowly rotating c l ouds. A s shown in lCha fc Whitworth I (120031 ) , iMatsumoto &: Hanawa 
(120031 ). and iMachida et al. I (l2005bl ). fragmentation occurs in a collapsing cloud core and 
binary or multiple stars are formed in initially rapidly rotating clouds. When binary or 
multiple stars are formed, the angular momentum is redistributed in the orbital and spin 



angular momenta of each star (iMachida et al. Il2005bl ). Thus, the binary or multiple star 



formation process is also important for the angular momentum problem. We will investigate 
the evolution in initially rapidly rotating clouds in subsequent papers. 



5.3. Comparison with Observations 



In our calculation, the formed protostars have magnetic fields of -Bprotostar = 0.16 kG 
(Model 16) to 1.53 kG (Model 32) in models with t] = 1 (non-ideal MHD Models 16, 22, 
24, 25, 26, 27, and 28). On the other hand, a surface magnetic field of the order of 1- 
3kG has been derived from Zeeman broadening measurements of CTTS photq spheric lines 
Jjohns-KruU et al. Ill999al lbl. I2OO1I : Icuenther et al"~lll999l : [j^vier et al. \\200^ . Thus, the 
magnetic fields derived from our calculation are consistent with observations. On the other 
hand, the rotation periods of protostars in our calculation are slightly shorter than the 
observations. The formed protostars have rotation periods of P = 0.15 days (Model 26) to 
1.53 days (Model 32) in models with Cr , = 1, while observed protostars have rotation periods 
of P = 1-10 days Jflerbst et al. Il2006h . 



Since our calculations show only the early star formation phase, in which the protostars 
have masses of ~ lO^^-lO^^M©. The magnetic field strength and angular momentum of a 
protostar may change in subsequen t gas accretion phase. T he magnetic field may be amplified 
by magneto-rotational instability (IBalbus fc Hawleylll99ll ) and by shearing motion between 
the star and the ambient medium. The angular momentum of a protostar may increase 
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due to accretion of gas with large angu lar momentum, and decrease by jets or star-disk 
interaction (e.g., Matt fc Pudritz 2005al lbl). To fully describe the magnetic field and angular 



momentum of protostars, we must calculate the cloud evolution up to the later accretion 
phase. However, our results indicate that at its formation epoch, a protostar already has 
a magnetic field strength and rotation period almost of the same order of magnitude as 
observations. Thus, the magnetic field and angular momentum of observed protostars may 
be largely determined by the early star formation phase before the main accretion phase. 

We have greatly benefited from discussions with T. Nakano, and K. Saigo. We also 
thank T. Hanawa for his contribution to the nested grid code. Numerical calculations were 
carried out with a Fujitsu VPP5000 at the Astronomical Data Analysis Center, National 
Astronomical Observatory of Japan. This work is supported by a Grant-in-Aid for the 
21st Century COE "Center for Diversity and Universality in Physics" from the Ministry 
of Education, Culture, Sports, Science and Technology (MEXT) of Japan, and partially 
supported by the Grants-in-Aid from MEXT (15740118, 16077202, 18740113,18740104). 
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Table 1: Model parameters and calculation results 
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Fig. 1. — Resistivity rj (left axis) and magnetic Reynold s number Rem ( right axis) as a 
function of the number density. The resistivity estimated by lNakano et al. I (120021 ) is plotted 
as a thick solid line (cj^ = 1) and the corresponding magnetic Reynolds number as a thin solid 
line. The magnetic Reynolds number is analytically acquired by using the free-fall velocity 
Vf and Jeans length Aj estimated for the given density and resistivity as Re^n = VfXji]'^. 
Resistivities multiplied by factors of 0.001 (c^ = 0.001), 0.01 (c^ = 0.01), 0.1 (c^ = 0.1), 
and 10 {cr^ = 10), and the corresponding magnetic Reynolds numbers are also plotted. The 
magnetic field is well coupled with the gas in the "coupled" region, while the magnetic field 
is decoupled from the gas in the "decoupled" region. Below the vertical broken line Re^ = 1, 
Ohmic dissipation is effective. Thus, in the shadowed region (3 x 10^^ ^ ^ 5 x 10^^ cm^'^), 
the magnetic field is effectively dissipated by Ohmic dissipation for = 1. 
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Fig. 2. — Evolution of the central magnetic field (upper panel) and angular velocity Qc 
(lower panel) against the number density at the center of the cloud for Models 3, 8, 16, 22, 
and 23. The rotation periods P — 27r/flc are also plotted on the right-hand axis of the lower 
panel. 
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1/2 

Fig. 3. — Magnetic field normalized by the square root of the density pc in the unit 
of initial sound speed Cg.o plotted against the number density at the center of the cloud for 
Models 1-6 (upper panel) and 7-11 (lower panel). 



Fig. 4. — Density (color-scale) and velocity distribution (arrows) on the cross-section in 
the z = plane (upper panels) and y = plane (lower panels) for Model 12 [{a,uj,Cr,) = 
(0.01, 0.01, 0)]. Panels (a) through (f) are snapshots at the stages (a) Uc = 1.2 x lO^cm"'^ 
{I = 4; initial state), (b) 1.0 x 10^ cm"^ (/ = 9-11; isothermal phase), (c) 3.4 x 10^^ cm~^ 
(/ = 15, 16; adiabatic phase), (d) 2.5 x lO^^cm"^ (/ = 17-20; second collapse phase), (e) 
1.8 X lO^^cm-^ (/ = 24, 25; protostellar phase), and (f) 4.4 x 10^^ cm'^ (/ = 24, 25; 
calculation end), where I denotes the level of the subgrid. The dotted lines indicate the first 
core (white) and second core (black), surrounded by the shock surfaces. The red thick hues 
indicate the border between the infalling and outflowing gas (fout > km s~^). The level 
of the subgrid is shown in the upper left corner of each upper panel. The elapsed time t, 
density at the center ric, and arrow scale are denoted in each panel. 
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Fig. 5. — Three-dimensional projections of Fig. HI The panels (a)-(f) have the same scale 
and epoch as panels (a)-(f) of Fig. HI The structure of the high-density region (p > O.lpc, 
red iso-density surface) and magnetic field lines (black-and-white streamlines) are plotted in 
each panel. The density contours (false color and contour lines), velocity vectors (arrows) 
on the cross-section in the x = 0, y = 0, and z =0 plane are, respectively, projected on the 
sidewalls of the graphs. The grid level /, elapsed time t, density at the center ric, and grid 
scale are denoted in each panel. 
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Fig. 6. — Evolution of models with different diffusivity rj. Upper panels (a) and (d): 

1/2 

Magnetic field B^c normalized by the square root of the density pc in the unit of initial 
sound speed Cgfi against number density at the center of the cloud for (a) Models 12-17 

and (d) Models 18-23. Middle panels (b) and (e): Plasma beta Sttc'Ip/B'^ within the region 
p > O.lpc for the same models as panels (a) and (d). Lower panels (c) and (f): J/M"^ in the 
unit of 47rG/cs,o within p > O.lpc for the same models as panels (a) and (d). 
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Fig. 7. — Same as Fig.Hlbut for Model 16. Panels (a) through (f) are snapshots at the stages 
(a) Tic = 4.9 X 10^3 cm-3 (/ = 15, 16), (b) 2.7 x lO^^cm-^ (/ = 17-20), (c) 3.4 x lO^^cm-^ 
(/ = 19-21), (d) 6.0 X lOi^cni'^ (/ = 20-23), (e) 8.1 x 10^^ cm-^ (/ = 21-25), and (f) 
1.3 X 10^°cm~^ (/ = 24, 25), where / denotes the level of the subgrid. The dotted lines 
indicate the first core (white) and second core (black), surrounded by the shock surfaces. 
The red thick lines indicate the border between the infalling and outflowing gas (fout > 
km s~^). The level of the subgrid is shown in the upper left corner of each upper panel. The 
elapsed time t, density at the center ric, and arrow scale are denoted in each panel. 




Fig. 8. — Same as Fig. [5] but for Model 16. Panels (a)-(f) are for the same scale and epoch 
as panels (a)-(f ) of Fig. [71 
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Fig. 9. — Velocity and plasma beta for Model 16 at the same scale and epoch as Fig. [7]g?. 
Upper panel {z > 0): Velocity of the 2;-component Vz (color-scale, contours) and velocity 
distribution (arrows) on the cross-section in the y = plane. Lower panel {z < 0): Plasma 
P Pp (color-scale, contours) and magnetic field (arrows) on the cross-section in the y = 
plane. 
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Fig. 10. — Evolutions of the magnetic field and angular velocity for various initial rotation 

1/2 

rates. Upper panel: Magnetic field Be normalized by the square root of the density pc in 
the unit of initial sound speed o against the number density at the centers of the clouds for 
Models 24-28. Lower panel: Angular velocity ilc normalized by the free-fall rate [{ATrGpcY^'^] 
at the centers of the clouds for Models 24-28. 
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Fig. 11. — (a) Spatial distribution of magnetic flux. The magnetic flux normalized by the 
critical mass $/$cri as a function of the cumulative mass from the center of the cloud at the 
initial state (broken hne) and at the end of the calculation for Models 29 (sohd hne) and 
30 (dotted line) are plotted, (b) Evolution of $/3>cri within p > O.lpc against the central 
number density for Models 29-32. (c) Spatial distribution of specific angular momentum. 
Specific angular momentum j as a function of the cumulative mass from the center of the 
cloud at the initial state (thick broken line) and the end of the calculation for Models 29 
(solid line), 30 (thin broken line), and 33 (dash-dotted line) are plotted, (d) Evolution of j 
within p > O.lpc against the number density at the center of the cloud for Models 29, 30, 
and 33. 



